Ultracold Rb-OH collisions and prospects for sympathetic cooling 



o 
o 



Manuel Lara,^ John L. Bohn/ Daniel Potter,^ Pavel Soldan,^ and Jeremy Hutson^ 

'JILA, NIST and Department of Physics, University of Colorado, Boulder, CO 80309-0440 
'^Department of Chemistry, University of Durham, 
South Road, Durham, DHl 3LE, United Kingdom 
^Doppler Institute, Department of Physics, Faculty of Nuclear Sciences and Physical Engineering, 
Czech Technical University, Bfehovd 7, 115 19 Praha 1, Czech Republic 
(Dated: February 2, 2008) 

We have computed ab initio cross sections for cold collisions of Rb atoms with OH radicals. We 
predict collision rate constants of order 10~^^ cm^/s at temperatures in the range 10-100 mK at 
which molecules have already been produced. However, we also find that in these collisions the 
molecules have a strong propensity for changing their internal state, which could make sympathetic 
cooling of OH in a Rb buffer gas problematic in magnetostatic or electrostatic traps. 



PACS numbers: 34.20.Mq, 34.50.-s, 33.80.Ps 



Oh' 



O 
Oh' 



> ■ 

oo ■ 

o : 

o • 

^ : 

o . 

^ : 

Oh. 



X 



A great way to make something cold is to place it in 
thermal contact with something else even colder. Thus, 
for many years, it has been possible to cool one species 
of ion sympathetically, by placing it in contact with an- 
other species that is being actively cooled P, 01 • More 
recently, ultracold neutral atoms have also been sympa- 
thetically cooled 0. At slightly higher temperatures, 
atoms and molecules have been sympathetically cooled in 
a helium buffer gas Q . Molecules are widely regarded as 
worth cooling to ultralow temperatures (typically around 
1/iK), where de Broglie wavelengths become large, and 
intermolecular interactions are dominated by long-range 
forces. Under such circumstances, the possibility of per- 
manent electric dipole moments is expected to lead to 
novelties in such areas as precision measurement, colli- 
sions and chemistry, quantum degenerate matter, and 
quantum information theory |5l|. 

A host of molecular cooling techniques have been pro- 
posed 0, but two that stand out are buffer-gas cooling 
(BGC), in which the cold He gas cools the molecules, 
and Stark deceleration (SD), in which polar molecules are 
slowed by carefully designed time-varying electric fields 
0. These techniques tend to produce "lukewarm" gases, 
with translational temperatures in the 10-100 mK range. 
They are therefore somewhat at a disadvantage with re- 
spect to the direct production of molecules by photoas- 
sociation from a gas of ultracold atoms 0|- Nevertheless, 
BGC and SD are extremely appealing for the far greater 
variety of molecular species that they can cool. 

In this Letter we are concerned with the possibility 
of turning a 10 mK gas of molecules into a 10 /iK gas, 
by sympathetically cooling with Rb. For this cooling 
to be effective, the cross sections must be favorable. 
Namely, elastic scattering should occur frequently (with 
large cross sections ctoi), to make the necessary thermal 
contact. On the other hand, inelastic collisions that 
change the state of the molecules should occur infre- 
quently (with small cross sections (Tinei)- This is because 
either magnetostatic or electrostatic trapping demands 



that the molecules remain in a well-defined internal state. 
Until now, cross sections for m-level-changing cold col- 
lisions of alkali atoms with molecules were wholly un- 
known. In this Letter we present the first such collision 
results, using a set of complete ab initio potential energy 
surfaces (PES's), and incorporating the hyperfine struc- 
ture of the collision partners. 

Our prototype system for this study consists of OH 
molecules in a bath of Rb atoms. This is a particularly 
timely example: the OH radical has been successfully 
slowed by Stark deceleration techniques in at least two 
laboratories, which can now produce ~ 10 mK packets 
of these molecules on demand @, At the same time, 
Rb is easily cooled and trapped in copious quantities at 
tens or hundreds of /iK, making it an ideal target for 
the molecules. In addition, the Rb-OH collision system 
is subject to a "harpooning" process, similar to that in 
Rb-NH , where during a collision the valence electron 
jumps from Rb to OH [l(]j . This process is without prece- 
dent in the cold collisions literature. Thus, even apart 
from sympathetic cooling, the Rb-OH scattering system 
is quite rich from a cold collisions perspective. 

We begin by describing the potential energy surfaces 
(PES's) of the system, as functions of (i?, 9), where R is 
the distance between the atom and diatom, and 9 is the 
angle that the Jacobi vector R makes with the OH axis. 
The OH monomer has a ^H ground state arising from 
a TT^ configuration, while the ground state of Rb is ^S. 
This produces ^H and ^H states for linear RbOH, which 
split into ^A', ^A" and ^A', ^A" surfaces at nonlinear 
geometries. We will refer to these as the covalent states. 
In addition, RbOH has an ion-pair state analogous to the 
ones previously found for RbNH |^ . The ion-pair thresh- 
old for Rb+ CS) + OH- (iE+) lies only 2.35 eV above 
the ground state at i? = oo and produces a (^^') 
state that cuts steeply down, crossing the covalent states 
near 6 A. This is an actual crossing at linear geometries, 
where the ion-pair and covalent states have different sym- 
metries, but there are avoided crossings between the two 
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states of ^A' symmetry at nonlinear geometries. There 
is therefore a conical intersection between the two ^A' 
surfaces at linear geometries. 

Potential energy surfaces for all 5 electronic states of 
RbOH were calculated by complete active space self- 
consistent field (CASSCF) calculations followed by mul- 
tireference configuration interaction (MRCI). A full de- 
scription will be given elsewhere jllj. All calculations 
used the MOLPRO package [13, with aug-cc-pVTZ ba- 
sis sets in uncontracted form for O and H and the 
ECP28MWB small-core quasirelativistic effective core 
potential for Rb, with the valence basis set from 
Ref. 0. The CASSCF calculation included aU configu- 
rations that can be formed from the lowest (10,3) orbitals 
of (a', a") symmetry, with the lowest (5,1) orbitals doubly 
occupied. Energies were calculated in Jacobi coordinates 
{R, 9) on a grid of 25 unequally spaced points in R from 
2 to 12 A and 11 Gauss-Lobatto quadrature points in 
9. The two ^A' surfaces were transformed to obtain two 
quasidiabatic diagonal surfaces and a coupling surface, 
using a mixing angle derived from matrix elements of 
the OH Lz operator. Finally, the quasidiabatic covalent 
energies were extrapolated to i? = cx) using Cg and Cy 
coefficients obtained from coupled cluster (CCSD) calcu- 
lations at R — 15, 25 and 100 A. The singlet and triplet 
sum and difference surfaces were expanded in normalized 
associated Legendre functions and the radial coefficients 
were interpolated using the reciprocal-power reproducing 
kernel Hilbert space (RP-RKHS) method [lllll|. 

The resulting diabatic covalent surfaces have wells 337 
and 511 cm~^ deep at linear Rb-OH geometries for the 
singlet and triplet surfaces respectively. For the A' sur- 
faces this linear well is the absolute minimum, but the A" 
surfaces have absolute minima at non-linear geometries 
{9 « 125°) with depths of 405 and 615 cm^^ respec- 
tively. The ion-pair state is very much deeper, with a 
well ~ 26000 cm~^ below the neutral threshold. 

To perform scattering calculations on these surfaces, 
we expand the wave function into an appropriate basis 
set. This set will consist of the quantum numbers of the 
atom and diatom in the separated limit. For the atom, 
electronic (sa) and nuclear (ia) spins are coupled to make 
total spin |/a"^/a) in the lab frame. For the diatom, the 
electronic wave function is specified in Hund's coupling 
case (a) by the total electronic angular momentum j, 
with projections m in the lab frame and uj in the frame 
rotating with the molecule. As usual, ui is separated into 
its orbital and spin components, lu = X + a. Linear com- 
binations \Xide) — (|Aw)-|-e| — A — a;))/-\/2 are constructed 
to produce states of good parity appropriate to the zero- 
electric-field case we consider here. Finally, j is coupled 
to the nuclear spin id of the hydrogen atom to yield total 
diatom spin Ifciruf^). The partial wave quantum num- 
bers \LMl) account for the relative orientation of the 
collision pair in the lab frame. Thus the basis states of 



our calculation are given by 

\sdXCje{jid)fdmfJ\{saia)famfJ\LML). (1) 

In this basis the total parity p — e(— 1)-'^*''+^ and the 
lab-projection of total angular momentum, rud + rua + 
Mi^, are conserved quantities. In zero field, the total an- 
gular momentum = f^+fa-l-L would also be conserved; 
however, we anticipate considering the action of a field 
on the collisions, which would mix different values. We 
therefore do not exploit this symmetry here. 

Cast in this basis, the Schrodinger equation for Rb- 
OH collisions takes the standard form of a set of close- 
coupled equations in the variable R. To solve these, we 
propagate the log-derivative matrix Y = ip~^dip/dR us- 
ing a variable-step version of Johnson's algorithm p^ . 
Because of the multiple PES's and the inclusion of hy- 
perfine degrees of freedom, the total number of scattering 
channels required to solve the complete problem includ- 
ing spin (in excess of 25,000) is prohibitively large. 

To overcome this obstacle, we exploit a kind of frame- 
transformation procedure [l9|. We identify a suitable 
radius R ~ Rq, and define R < Rq and R > Rq a,s 
the "short-range" and "long-range" parts of the calcula- 
tion, respectively. For Rb-OH, we choose Rq = 17 ag. 
In the short-range region, the PES's are deep, and the 
hyperfine effects are small. We therefore compute Y in 
a basis of decoupled nuclear spins and neglect couplings 
between blocks of the Hamiltonian with different nuclear 
spin projections rriia or mid- In this pilot study, we have 
neglected the influence of the ion-pair channel. 

At long range, R > Rq, the full hyperfine structure is 
restored, and Y is transformed into the basis (1) for fur- 
ther propagation and matching to spherical Bessel func- 
tions to yield scattering matrices. In this region, collision 
channels in which the molecule is excited into higher- 
lying rotational or spin-orbit states are already strongly 
closed. We therefore eliminate these channels. In the 
inner region, we use partial waves up to L = 28 and ro- 
tational states up to j — 11/2. In the outer region, we 
employ partial waves up to L = 22, but with full hy- 
perfine structure of the j = 3/2 rotational ground state. 
Because of these approximations, no single calculation 
requires more than 2000 channels ■ Numerical checks 
suggest that the magnitudes of cross sections computed 
in this way are accurate to within a factor of ^ 2. 

For concreteness, we calculate cross sections for a beam 
experiment in which the incident wave vector is coin- 
cident with the laboratory quantization axis. Energy- 
dependent cross sections for two different initial states of 
the collision partners are shown in Fig. In the first 
case. Fig. Q^a), the Rb atom is initially in its "stretched" 
state, with |/a"^/„) = |22), and the molecule is in its 
higher-lying, / parity state, and also spin stretched, with 
l/d"^/d) = 1 22). Both these states are weak- magnetic- 
field seeking, meaning that they can be trapped magnet- 
ically. In addition, OH in this state can be trapped elec- 
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FIG. 1: Rb-OH collision cross sections versus energy. Solid 
(dashed) lines represent elastic (inelastic) cross sections. In 
(a) is shown results for the incident channel fli3/2,fd'nifd = 
22,/)oH Ifarrifa ~ 22)Rb, where both the atom and diatom 
are weak-magnetic-field seeking, and the OH is also weak- 
electric-field seeking. In (b) is shown the incident channel 
\'^^3/2, fdrufd = ll,e)oH Ifanifa = 1 - l)Rb. Again, both 
species are in weak-magnetic-field seeking states, but the OH 
is now in a strong-electric-field seeking state. The points in- 
dicate the Langevin cross section, Eqn. (2) 



trostatically. Solid lines denote elastic cross sections CToi, 
in which all internal quantum numbers are retained after 
the collision, and dashed lines refer to inelastic cross sec- 
tions (Tinci, which represents the sum of partial cross sec- 
tions to many possible outcomes distinct from the initial 
channel (including nia and m^- -changing collisions that 
do not release energy at zero field). 

For collisions between S-state atoms, being in a 
stretched state greatly reduces inelastic cross sections. 
This is simply because the m/^ projection of one atom 
cannot be reduced at the expense of raising rrif^ of the 



other atom. Rather, to conserve angular momentum, 
such a collision can occur only if the partial wave projec- 
tion Ml can change, which requires an anisotropic inter- 
action. For two S-state atoms, the only anisotropic inter- 
action is the magnetic dipole interaction between valence 
electrons, which is very weak. For Rb-OH collisions, by 
contrast, the PES itself is highly anisotropic, providing a 
ready avenue for spin-changing collisions. Indeed, inelas- 
tic scattering is as likely as, or more likely than, elastic 
scattering over much of the energy range. This trend is 
especially pronounced at low energies (below ~ 1 mK), 
where the Wigner threshold law sets in. The elastic cross 
section must tend to a constant at threshold, but ex- 
oergic inelastic processes must yield cross sections that 
diverge as 1/Ve. Note that the actual magnitudes of 
cross sections in this limit are strongly subject to details 
of the PES, details which are typically uncovered only 
by experiments. Thus the cross sections in Fig. ^ while 
suggestive, cannot be taken as definitive predictions. 

At higher energies, many partial waves can contribute 
to the scattering, and the threshold laws are less signifi- 
cant. Because the electron spin of the OH is strongly tied 
to the rotating molecule, spin orientation in the lab frame 
is completely disrupted during a collision. This process 
can be approximated simply by a Langevin model |lfll |: 
for a given energy i?, there is a maximum partial wave 
L{E) for which E lies above the centrifugal barrier of the 
long-range potential -Ce/R^ + ffL{L + \)/2^R^. Here 
Cg = 325£^fran is the isotropic Van der Waals coefficient 
of the PES [ll|, and ^ the reduced mass of the colli- 
sion pair. In the Langevin model, we assume that any 
partial wave smaller than L{E) contributes to inelastic 
scattering with unit probability, i.e., any time the atom 
and molecule are near one another, they are essentially 
guaranteed to disrupt the electron's spin. This idea leads 
to a cross section 
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This cross section is also shown in Fig. ^a) , by the solid 
symbols. It is clearly getting the trend and the order of 
magnitude of the cross sections correct. 

We therefore predict that Rb-OH elastic cross sections, 
at temperatures of tens of mK typical of Stark deceler- 
ators, will be on the order of 10~^^ cm^. At the cor- 
responding collision velocities of several hundred cm/s, 
this corresponds to elastic collision rate constants on the 
order of 10^^^ cm'^/s, which is easily large enough to 
provide sufficient thermal contact to allow sympathetic 
cooling in overlapped traps of OH and Rb. However, 
inelastic rate constants are of the same order of magni- 
tude, which bodes ill for sympathetic cooling, at least in 
magnetostatic or electrostatic traps. Including the ion- 
pair channels would likely only disrupt the internal states 
more severely, and would not change this basic conclusion 
about sympathetic cooling. 
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We turn now to an alternative pair of initial states, 
whose scattering cross sections are shown in Fig. ^ (b) . 
Here the atom is in its lowest magnetically trappable 
state, with |/a™/„) — |1 ~ l)j ^^nd the molecule is in 
its lower-energy e state, with total spin Ifdm/a) — 
(also magnetically trappable). At zero field, both CTci and 
(Tinei are energy-independent at low energies, since all exit 
channels are iso-energetic with the incident channel. In 
nonzero field, the |11)oh|1 — l)Rb channel lies above other 
m-components in energy, and CTinoi would again diverge. 

By the same token, sympathetic cooling is possible at 
low energies for OH molecules in the |1 — 1) state and Rb 
atoms in the |11) state. This is the lowest-energy chan- 
nel of all: there are no inelastic channels energetically 
available and (Tinoi will vanish. These states are not mag- 
netically trappable, but could be confined in an optical 
or microwave 20] dipole trap or an alternating current 
trap [iJl . A small static magnetic field will be needed to 
maintain the projection quantum numbers. 

At higher collision energies, other forms of inelastic 
scattering become energetically allowed. Once the colli- 
sion energy surpasses the height of the p-wave centrifugal 
barrier (at about 1.6 mK), the molecule can shed angu- 
lar momentum into the partial wave degree of freedom. 
At this point (Tinei climbs to a large value, comparable to 
(Tel- Further, at about 4mK, the / = 2 hyperfine state 
of OH becomes energetically allowed, providing another 
route to inelastic collisions. By examining the partial 
cross sections, we find that these two avenues for inelas- 
tic collisions are roughly equally likely. At still higher 
collision energies, (Tinci is again roughly approximated by 
the Langevin result Q , although there is now more struc- 
ture, owing to a large number of Feshbach resonances to 
fine-and hyperfine excited states. 

These results present a cautionary tale for sympathetic 
cooling using alkali atoms. One may regard inelastic 
scattering as a "Murphy's Law" process |22|. Namely, 
if the internal state of the molecule can change to pro- 
duce an unfavorable result, it will do so. The key to 
making sympathetic cooling viable lies in eliminating un- 
desirable channels as far as possible. For example, a light 
collision partner with a small Ce coefficient will produce 
a centrifugal barrier at higher energy, preventing partial 
waves from accepting angular momentum. This circum- 
stance explains, at least partly, the ability of a He buffer 
gas to cool molecules without badly disrupting spin ori- 
entation ^23,, ,24, In the present context of alkali 
atoms, consider Li, which is about 12 times lighter than 
Rb, and about half as polarizablc. Its p-wave centrifu- 
gal barrier upon colliding with OH is on the order of 10 
mK. In addition, a molecule that is better described by 
Hund's coupling case (b), in which the electron's spin is 
only weakly coupled to the molecular rotation axis, may 
help weaken CTinci below the Langevin limit [25l |. 

In summary, we have performed the first ah initio 
scattering calculations for an open-shell, ground state 



molecule colliding with an alkali atom, incorporating the 
hyperfine structure of both collision partners. The results 
suggest that elastic cross sections are sufficiently large for 
sympathetic cooling to occur, yet equally large inelastic 
cross sections probably hinder this application. Nonethe- 
less, the results do not preclude the possibility that, by 
applying electric or magnetic fields, inelastic cross sec- 
tions could be suppressed [2^|23|. This will be the sub- 
ject of a future study. 
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